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The evolution of two-component cold atoms on a ring with spin-orbit coupling has been studied 
analytically for the case with N noninteracting particles. Then, the effect of interaction is evaluated 
numerically via a two-body system. Two cases are considered: (i) Starting from a ground state the 
evolution is induced by a sudden change of the laser field, and (ii) Starting from a superposition 
state. Oscillating persistent spin-currents have been found. A set of formulae have been derived to 
describe the period and amplitude of the oscillation. Based on these formulae the oscillation can be 
well controlled via adjusting the parameters of the laser beams. In particular, it is predicted that 
movable stripes might emerge on the ring. 
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I. INTRODUCTION 

It is well known that the study of the motion of charged 
particles under a magnetic field is an essential topic in 
both macroscopic and microscopic physics. In particular, 
a number of distinguished quantum mechanic phenom¬ 
ena, such as the Aharonov-Bohm (A-B) oscillation and 
the fractional quantum Hall effect (FQHE), are caused 
by the magnetic gauge field. [H-I^ After the experimen¬ 
tal realization of the condensation of neutral atoms with 
nonzero spin,Q a great interest is to create a light- 
induced gauge vector field to achieve the spin-orbit cou¬ 
pling (SOC) so that various magnetic-electronic phenom¬ 
ena of charged particles occurring in condense matters 
can be ciroied into the world of condensates of neutral 
atoms.Since 2005 there are proposals by the the¬ 
orists for producing the light-induced gauge vector field 
for the condensates with multi-component neutral cold 
atoms. (Tol - fTsjl The first experimental realization of the 
SOC in condensates was implemented in 2009 by dressing 
the cold atoms with two counter propagating polarized 
laser beams. [il[i3 This technique opens a new perspec¬ 
tive in the field of BEG, more and more experimental [iSl- 
[2l| and theoretical [2a - [^ results are reported. It was 
found that the total magnetization (or spin-polarizability, 
which measures the difference of the densities of the spin- 
components) varies with time. It implies the existence 
of spin-currents. [H, m, Besides, the collective-mode 
dynamics has been studied, the spatial motions along 
different directions could be correlated (say, a dipole 
mode might induce a breathing mode in its perpendic¬ 
ular plane). The instability caused by the collective 
mode has also been studied.!^ 
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On the other hand it is now possible to trap a con¬ 
densate in ring geometry. Long-lived rotational super- 
flows have been induced in this kind of systems. [20l l25l - 
The creation of the spin-currents of two-component 
condensates on a ring, due to the introduction of the laser 
beams, has been experimentally realized. It was found 
that the stability depends strongly on the initial ratio of 
the two components, However, the oscillation of the 
spin-currents on the ring has not yet been observed. 

Obviously, the study of the cold spin-currents (i.e., the 
current of each spin-component) caused by the SOC is 
just beginning. It is expected that the oscillation of 
the spin-currents might also emerge in the ring geom¬ 
etry (because a circular motion is equivalent to a linear 
motion with a periodic boundary condition). To con¬ 
firm this suggestion, a one-dimensional model of the two- 
component condensates on a ring under SOC is adopted 
in this paper. The aim is to clarify from the theoretical 
aspect how the details of oscillation (the period and am¬ 
plitude) depend on the parameters of the laser beams. 
We believe that the interaction might not be crucial. 
Therefore, the interaction is neglected firstly so as to 
emphasize the decisive role of the laser beams. Then, 
the effect of the interaction is evaluated via a two-body 
system. 


II. TIME-DEPENDENT SOLUTION 

Let the quasi-spin s be introduced to describe the two 
components of an atom as usual. The atomic state with 
Sz = 1/2 (—1/2) is named the up- (down-) state. When 
two counter-propagating and polarized laser beams are 
applied, due to the SOC, the two components of atoms 
might move towards opposite directions along the beams 
and can transform to each other via a spin-flip. When 
the interaction is neglected, the Hamiltonian is just H = 



2 


Ei/ii, where hi is for the i-th particle. Neglecting the 
subscript i, 


h = 




( 1 ) 


where 0 is the azimuthal angle along the ring.Q The 
unit of energy is hereafter i?unit = h‘^/{ 2 niR^), where 
m is the mass of an atom, R is the radius of the ring. 
/3 = koR, where 2 fco is the momentum transfer caused 
by the two lasers. 7 = , where r 2/2 is the strength 

of Raman coupling causing the spin-flips accompanied 
by the momentum transfer. Let £spiit be the Zeeman 
energy difference between the two spin-states, and ujg is 
the frequency difference between the two laser beams. We 
define the Raman detuning 5 = (ffspiit — hujs)/i‘2Eu^ait)- 
In addition to 7 , 6 is an important quantity because it 
can be tuned so that the related process could be energy- 
adapted (see below). Note that, due to the ring geometry, 
2/3 must be an integer. It implies that the transfer of 
momentum will be suppressed unless the transfer is close 
to a specific set of values depending on R. Thus, an 
additional constraint is required under the ring geometry. 

To illustrate the connection with conventional SO cou¬ 
pling, we define a t/-transformation so that the up-state 
(down-state) is multiplied by a factor (e“*^®). Then, 
by introducing the Pauli matrices Uz and ov, UhU~^ can 
be written in a more familiar form as [§, [ij, EBi 113) Hsjl 


UhU ^ + (2) 

do 

where ai is just a unit matrix with rank 2 . 

h has two groups of eigenstates, they can be written 

as [HillHI 

= s\n{pk)‘^k-0,^ + cos{pk)(pk+0,i, (3) 

= cos(pfe)(pfc_/3,t - sin(pfe)(^/c+/34, (4) 

where 

k is an integer (half-integer) when /3 is an integer (half¬ 
integer), sinpfc = s^v^(afe^^^fc)3T"d)7(2^, cos pk = 
{ak + 2kp — <5)/(2afe), Sry is the sign of 7 , ak = 
Vi^kp - (5)2 -I- 7 ^, and pk is ranged from —tt/2 —^ 7 r/ 2 . 
One can see that the eigenstates have two notable fea¬ 
tures: 

(i) Each of them contains both the up and down com¬ 
ponents, they are propagating with angular momenta 
k — P and k + P, respectively. Where 2/3 is propor¬ 
tional to the momentum transfer caused by the counter- 
propagating laser beams. 


(ii) In each the weights of the two components 
can be controlled by tuning S. When 6 < 2kp, we have 
— 7 r /4 < Pk < 7 r/ 4 , and therefore the weight of the up- 
component is larger in ^ but smaller in . Whereas 
when 6 > 2 kp, 0 ^ ^ will contain more down-component. 
In particular, when S is tuned so that 6 = 2kp, the 
weights of the two components in all the eigenstates are 
equal. 

The eigenenergies of 0^^^ are = k'^ + P'^ ± ak- 
Obviously, eIt'’ < 

The time-dependent solution ijj{9,t) of the single¬ 
particle Schrodinger equation starting from an initial 
state 0 init can be formally written as 

0(6»,t) = e“*^''0i„it 

= (7) 

k\ 


where r = tEunit/h (for ®^Rb and R = 12p,m as given 
in [ 2 ^, t = 0.398rsec), A = ±, fc runs over the (half-) 
integers when /3 is an (half-) integer. 


III. EVOLUTION CAUSED BY A SUDDEN 
CHANGE OF THE LASER FIELD 


In the experiment reported in (2lj | a strong oscillation 
of the magnetization induced by a sudden change of the 
laser field was reported. It is assumed that a set of pa¬ 
rameters /3, 71 , and i5i (the first set) are given initially, 
and the initial state is just the ground state (g.s.) of h as 

0init = 0gs = cos(/9^^G35g-/3,t “ where q 

depends on the three parameters so that the associated 
energy Eq~'’ is the lowest. The superscript in implies 
that it is obtained from the first set of parameters. Then, 
7 i and di are suddenly changed to 72 and 82 (the second 
set). Accordingly, we have a new Hamiltonian and a new 
set of eigenstates. With them ip{9,t) becomes 

0(0, f) = - fd^q+p,i), (8) 

where /„ = cos[p^P) cos{a^q\) icos{ 2 p^q^ — 

p^q^)s\a{af\), and fa = s\-o{p^P) cos[a^q\) + 
jsin(2pg^^ — sin(ag^V). Where a^q'’ and pg^^ 
are obtained from 72 and 82 - From 0(0,t), we know 
that the two components of the g.s., (pq-p^-f and 
<Pg_P 7 , remain unchanged. However, the coefficients of 
composition are no more constants but oscillate with r. 

Let the time-dependent densities of the up- and down- 
component be defined from the identity 0 “''( 0 , t) 0 ( 0 , t) = 
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n-^( 9 ,T) + ni( 0 ,T). Then, we have 

n-f = ^{cos^(p^^)) + [cos^(2p(2)-pW)_cog2(^(i))] 

X sin^(a(^V)}, (9) 

rii = ^{sin2(p«) + [sin2(2pf-p«)-sin2(pW)] 

X sin^(a^^V)}. (10) 

The magnetization (or spin-polarization) is defined as 

Pz = (n-f - ni)/(n-i + ni). We have 

Pz = cos(2pW) + [cos( 4 p( 2 ) _ 2pW) - cos(2pW)] 

X sin^(a^^V). (11) 

This formula gives a clear picture of a harmonic 9 - 

/o\ 

independent oscillation with a period Tp = nfaq and 

an amplitude Aamp = | cos( 4 pq^^ — 2p[^'^) — cos(2pq^^)|. 

Based on the second set of parameters let us de¬ 
fine the energy difference of the two components as 

^diff = - {‘Pg+P,i\h^'^M‘Pg+l^,i)\ = 

{ 2 q$ — (52). Then, the frequency of oscillation l/r^ = 

i + 72- This formula demonstrates that both 

the strength of the Raman coupling and the energy dif¬ 
ference are crucial. A larger I72I leads always to a higher 
frequency > |72|/7r. This frequency can be further tuned 
by varying ^2 around 2 q/ 3 . In particular, when ^2 is tuned 
so that is zero, l/rp will arrive at its conditional 

minimum |72|/7r. Since depends on q but not on 
the other details of the initial state, l/r^ depends essen¬ 
tially on the second set of parameters when q is given. 

On the other hand, the amplitude Aamp depends on 
both sets of parameters explicitly. When the first set of 
parameter has been given so that is fixed, Aamp will 
arrive at its conditional maximum at a specific value of 
namely, p^^'* = p''q'‘ /2 (if 72 < 0 and p[^^ < -7r/4), 
or = p^q'’ /2 + 7r/2 (if 72 > 0 and p^q^ < —7r/4), or 
= p^q^ 12 — 7r/4 (if 72 < 0 and —7r/4 < p^q^ < 7 r/4), 
or = Pq V2 + Tr/d (if 72 > 0 and —7 r/4 < p^q^ < 7 r/4), or 
= p^q^ 12 — 7 r/2 (if 72 < 0 and p^q^ > 7 r/4), or = 

(if 72 > 0 and p^q^ > 7r/4). It implies that, no matter 
how the hrst set of parameters are, one can tune the sec¬ 
ond set so that the amplitude is conditionally maximized, 
namely, Aamp = 1 + | cos( 2 pg^^)|. In particular, when 
7i = 0, we have p'q^ = ±7r/2 or 0. In this case, if 72 0 

/o') 

and 61 is tuned so that = 0, then p\ ’ = (72/|72|)7’‘/4 
and the amplitude would arrive at its absolute maximum, 
namely, Aamp = 2 . Note that the condition p^q^ = ±7r/2 
or 0 implies that the initial state is pure without mix¬ 
ing (namely, either V'init = or P5-1-P.4.), while the 

fact Aamp = 2 implies a complete transformation of the 
two components (namely, from a pure up-state to a pure 
down-state, or vice versa). Therefore, 71 = 0 is required 
so that the maximal amplitude could be realized. On the 
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FIG. 1 : The frequency 1 /rp versus and 72. The oscillation 
is caused by a sudden change in 7 and < 5 . /3 = 3 and q = — 3 
are given. The values of the contours are in an arithmetic 
series. The contour closest to the right side has the largest 
value l/vp = 14 . 5 . 



other hand, in any cases, when 72 and 71 have the same 
sign and p^q^ —>■ pg^\ or 72 and 71 have opposite signs and 
p^q^ —7> p^q^ ± 7 r/ 2 , Aamp “0 and the oscillation damps. 
This is obvious from the expression of Aamp- 

In addition to the densities, one can further de¬ 
fine the current based on the equation of continuity as 
— ^(' 0 +(d,t)'!/)( 0 ,t)) = From this definition we 

found that the current j = j-|- -I- jj, in which = 
junitti-|-(q — /3) is the current of the up-component and is 
named up-current, while jj, = junitni{q + 13 ) is the down- 
current, where the unit is junit = The cur- 

rents are also oscillating with the same period ir/aij , and 
they are also 6>-independent. Obviously, when \j 3 \ > |q|, 
and jj, will have different signs. Therefore, counter 
propagating currents emerge. This is a distinguished fea¬ 
ture. 

The above formulae arise from a single-particle Hamil¬ 
tonian. For N-particle systems, when all the particles 
fall into the same state initially and the interaction 
is neglected, it is straight forward to prove that the above 
formulae hold also. 

To give numerical results, the radius is given at i? = 
12 pm in this paper. A contour diagram for the frequency 
1 /Tp versus 72 and <52 is shown in Fig.I. Where the 
minimum is located at 72 = 0 and E)^^g = 0 (namely, 
62 = 2 qf 3 = — 18 ), at which the oscillation vanishes. 

Two contour diagrams for Aamp versus 72 and <52 
with 7i and being fixed are shown in Fig. 2 , where 

all the contours appear as straight lines. This arises 
(2) 

from the fact that pq ' can be rewritten as a function of 

\l2/E)^^g\ = |72/(^2 — 2 qf 3 )\ together with the two signs 

( 2 ) 

of 72 and ^2 — 2 qf 3 (it implies that a given pq ' is associ- 
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ated with a straight line on the 82 — 72 plane). In Fig.2a 
(b) the given values of /3, 71 , and lead to g = —3 (3) 
and = 0 . 4357 r (O.OGSStt). According to the discus¬ 
sion given previously and when 72 > 0 , the maximum 
of Aamp should appear when = p^q'^ 12 = 0 . 2187 r in 

Fig.2a, and p^q^ = p^q^ I2 -\-tt j A = O.SIStt in 2b. This max¬ 
imum is associated with the straight dotted line marked 
in the figure (Incidentally, ii p\ — 7 r/ 4 , then the associ¬ 
ated straight line would be vertical. Since the two dotted 
lines have their p^q'^ close to 7 r/ 4 , they are nearly verti¬ 
cal). Whereas when p^q'^ = the associated straight 
line is marked by a solid line in the figure (where the 
small circle marks the place that 72 = 71 , and 82 = ^i). 
Three notable features are reminded: 

(i) The frequency does not depend on the first set of 
parameters, except q. A larger I 72 I leads always to a 
higher frequency which can be further tuned by varying 
82 around 2g/3. When <52 = 2qj3, the frequency arrives at 
its conditional minimum |72|/7r. 

(ii) Aamp depends on the ratios 71/^1 and 72 /^ 2 - When 
72 and 82 are given along the solid line (the slope of this 
line depends on 71 /di), the sudden change of 71 —>■ 72 
and (5i —^2 can not cause an evolution (i.e., Aamp = 0 ). 
When 72 and ^2 are given along the dotted line (also 
depends on 71 /^ 1 ), the amplitude will remain unchanged 
and is conditionally maximized, namely, Aamp = I -b 
I cos( 2 / 9 q^^)|. In any cases, one can always tune 82 so 
that the amplitude is either suppressed or conditionally 
maximized. 

(iii) Due to the fact that the contours are straight lines, 
when 72 is small (^ 0 ) there is a narrow domain of <52 
surrounding 2g/3 in which Aamp is highly sensitive to <52 
and appears as a sharp peak. At the peak the amplitude 
is large while the associated frequency is low. 

IV. EVOLUTION INITIATED FROM A 
SUPERPOSITION STATE 

Recently, the initial state was prepared in a superpo¬ 
sition state in an experiment of cold atoms with the ring 
geometry. In this experiment an additional radio fre¬ 
quency field was used for preparing the initial state 

•^init = sm{(t)/ 2 )ipq^ + cos{ 4 >/ 2 )ipq^, ( 12 ) 

where (j) determines the initial ratio of the two compo¬ 
nents and is tunable (0 < <() < tt). g is a tunable integer 
and marks the initial angular momentum of the particle. 
In this case we have the set of parameters 4>, g, /3, 7 , and 
8 . The corresponding time-dependent solution is 

+ fdi^q+2l3,i + fd2^q.\)^ ( 13 ) 

where = sin( 0 / 2 )e“^*'J^'^[cos(aq+/ 3 T) -b 

i cos{ 2 pq+p) sin(a,+/ 3 T)] 

! fu2 — 



FIG. 2: The amplitude Aamp versus 82 and 72 . I3 = 3, 
71 = 10, and 81 = 5(—5), are given in a (b). Accord¬ 
ingly, q = —3 (3) in a (b). The oscillation is caused by 
a sudden change from 71 to 72 and to ^ 2 . The dot¬ 
ted line marks the locations where = P^<p 12 (a), or 

p)f^ = p^q^ 12 -b 7 r /4 (b). Thus this line marks the condi¬ 
tional maximum of Aamp = 1 -b | cos(2pq^^)| = 1.917 in both 
(a) and (b), which is the largest value in each panel. The 
contours going away from this line are in an arithmetic se¬ 
ries and decrease to zero. The black small circle marks the 
location where 82 = (5i and 72 = 71 - The solid line passing 
through the circle marks the locations where p^^^ = pq^'^ and 
accordingly Aamp = 0 . 


-icos((/)/2)e2*«^^sm(2p,_/3)sm(a5_/3T) , fd^ = 

-isin((()/ 2 )e“^*'^^’‘sin( 2 pg+; 3 )sin(ag+; 3 r), and = 

cos((/)/2)e^*^^'^[cos(ag_^r) — icos{2pq-p)sm{aq-pT)]. 

Accordingly, we can obtain the up- and down-densities 
n-|-(0,T) and n-|-(0,T) given in the appendix. From them 
the magnetization can be obtained. It is apparent 
that, when (j) = tt or 0 (i.e., ' 0 init is a pure up-state or 
down-state) the magnetization has a very simple form as 

Pz = cos{Apq+p) -b (I - cos{Apq+p)) cos^iaq+pr), (14) 
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if (/) = TT, or 

Pz = - cos( 4 p,_/ 3 ) - (1 - cos(4p,_^)) cos^{aq-pT), (15) 
ii (j) = Q. 

These formulae also give a clear picture of harmonic 
oscillation, but the periods are different for the two cases 
of (j). 

When (j) = TT, Tp = Trjaq+p. If 6 is tuned so that S = 
2{q + /3)P = 6o, the frequency will be minimized, namely, 
(1/Tp)min = Wl\- When 5 goes away from So, the fre¬ 
quency increases. The amplitude Aamp = 1 —cos(4pq+^). 
Obviously, when 5 = So, we have Pq+p = ± 7 r/ 4 , and 
the amplitude would arrive at its maximum, namely, 
(^amp)max = 2 . Further more, when 7 — >■ ±0, one can 
prove that Pq+p will tend to 0 or ± 7 r/ 2 . In any of these 
cases Aamp = 0 and the oscillation disappear. This im¬ 
plies that 7 is the source of the oscillation. 

When 4> = Q, the discussion in the preceding paragraph 
holds also, except that q + P should be changed to q — 13, 
and So should be replaced by 2{q — (3)[3. An example 
of Aamp is shown in Fig.3. Fig.3 is of a similar form 
to Fig.2 (all the contours are straight lines converge to 
a point with 7 = 0 ) except that the dotted line in the 
former is exactly vertical. 

When (j) is neither tt nor 0, the oscillation is no more 
harmonic. In particular, both the up- and down-densities 
depend on 0 (refer to the appendix). It implies that the 
densities are not uniform on the ring as before. Conse¬ 
quently, stripes emerge along the ring, and the stripes 
move with time. An example of the evolution of Pz is 
shown in Fig.4. Where the relation S = 2{q + P)/3 holds. 
Accordingly, the curve with (j) = tt has the maximized 
Aamp = 2 and a minimized frequency. Whereas the curve 
with (f> = 0 has a smaller Aamp and a larger frequency. 
For (p = 7 r/ 2 , the anharmonic oscillation is clearly shown. 

An example of the up-density versus 9 with p = n/2 is 
shown in Fig.5, where the stripes emerge along the ring. 
Note that two waves 1^,7 and ifq- 2 p,^ are contained in 
the up-component, and ipq^i and (pq+ 2 p,i in the down- 
component. The stripes arise from the interference of 
these waves. Note that one of these waves contains an 
additional factor Accordingly, the stripes appear 

periodically with 9 as shown in Fig.5 and the number 
of peaks (valleys) in the domain (0, 27r) is 2/3 (say, the 
number is four in Fig.5, where (3 = 2). The locations of 
the peaks of the stripes move with time, the magnitudes 
of the peaks also change with time. 

As before, when all the particles are given in the same 
ipi-ait initially and the interaction is neglected, the above 
formulae and discussion hold also for N-particle systems. 

V. EVALUATION OF THE EFFECT OF 
INTERACTION 

When the interaction is taken into account, the to¬ 
tal Hamiltonian is H = where Vij = 

gS(9i — 9j), g = g^ (if both atoms are up), (both 



FIG. 3 : The amplitude Aamp versus 5 and 7 for the case that 
the initial state is a superposition state, g = 3 and {3 = 2 are 
assumed, p = 7r(0) are given in a (b). The dotted line marks 
the maximum of Aamp = 2 . The locations with 7 = 0 have 

Aamp — 0. 


down), or g^^ (one up and one down). In order to eval¬ 
uate the effect of interaction in the simplest way, we 
study a two-body system. Firstly, a set of single par¬ 
ticle states ipkfi = where = (J) or (°), 

and — fcmax ^ k < /cmax, are adopted. Accordingly, we 
have 2(2fcmax + 1) single particle states, and they are 
renamed as Based on ipj a set of basis func¬ 

tions for the two-body system = S[ipj{l)ipji(2)\ are 
defined, where S is for the symmetrization and normal¬ 
ization, and j > j'. 

Note that the strength of a pair of realistic Rb atoms 
is ffRb = 7.79 X lt)~^^HzC'rn?' (the differences in the 
strengths between the up-up, up-down, and down-down 
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FIG. 4: The magnetization Pz versus r observed at 0 = 0 for 
the case that the initial state is a superposition state with a 
0 marked by the curve, g = 3, /I = 2, 7 = 15, and 5 = 20 
are assumed. When (^ = tt or 0 the oscillation is harmonic. 
Otherwise, it is not. 



FIG. 5: Up-density n-i- versus 6 with 0 = 7 r/ 2 . The densities 
are given by 20 curves for a sequence of r. Starting from 
T = 0, in each step r is increased by 1/40. The 20 curves in 
the sequence are divided into 4 groups. The n-f of each group 
has been shifted up by 2 relative to its preceding group, for 
a clearer perspective. The five n-f in each group are plotted 
in solid, dash, dash-dot, dash-dot-dot, and dot lines, respec¬ 
tively, according to the time-sequence (r = 0 is the horizontal 
solid line in the lowest group). The parameters are q = 1, 
/? = 2, 7 = 10, and 5 = 1. The vertical dotted lines are also 
for guiding the eyes. 


pairs are very small and are therefore neglected). Since 
the atoms in related experiments are not distributed ex¬ 
actly on a one-dimensional ring but in a domain sur¬ 
rounding the ring, the effect of the diffused distribution 
should be considered. Hence, for each we define its 3- 

dimensional counterpart 

where the Gaussian function describes the dif¬ 



FIG. 6 : Pz of a 2-body system of Rb atoms versus r with 
interaction Vij = g5{9i — 9j), g = 0 (solid), lOOgefi (dash), 
and 2000(/eff (dash-dot). The initial state has 0 = tt (a) and 
= 0 (b). r-in = R/vTs. The other parameters are the same as 
in Fig.4 


fused distribution and measures the width of the dis¬ 
tribution. With we define an effective strength g^s 
so that for any pair of matrix elements 


ffeff J 9j)ipk3ii3fkitj.4 



Then, we have ges = ffRb which is the effective 

strength adopted in our calculation. 

It is assumed that both atoms are given at the same 
superposition state with </> = tt or 0 initially. When the 
Hamiltonian is diagonalized in the space expanded by , 
the eigenenergies Ei and eigenstates 'i’l = can 

be obtained, and the time-dependent state is 

^(0,t) (17) 

i 

where fPinit = '0init(l)^init(2). From '^/{9,t) the densities 
and Pz(t) can be calculated. An examples with fcmax = 
14 is shown in Fig .6 where the strength g is given at three 
values. When fcmax is changed from 14 to 12, there is no 
explicit changes in the pattern. It implies that the choice 
fcmax = 14 is sufficient in qualitative sense. 

In Fig .6 the curves ”1” and ”2” overlap approximately. 
It implies that the effect of interaction with a g < lOO^eff 
is very small. However, when g = 20005eff, the amplitude 
decreases with time explicitly as shown by the dash-dot 
curve, while the period remains nearly unchanged. Note 
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that, based on the Gross-Pitaevskii equation, the effect 
of the combined interaction imposing on a single particle 
from the other — 1 particles is similar to an appro¬ 
priate strengthening of the strength of the atom-atom 
interaction. Thus the explicit damp that happens when 
g = 20005 eff implies that, for an N-body system with a 
larger TV, the damp would actually occur. This is a topic 
to be clarified further. In fact, the damp of the oscil¬ 
lation has already been observed in all related existing 
experiments. [1^ 

VI. SUMMARY 

In summary, the oscillation of the cold atoms under the 
SOC and constrained on a ring has been studied analyti¬ 
cally without taking the interaction into account. Then, 
the effect of the interaction is evaluated numerically via a 
two-body system. Two cases, namely, the evolution that 
starts from a g.s. and is induced by a sudden change of 
the laser field, and the evolution that starts from a su¬ 
perposition state, are involved. The emphasis is placed 
on clarifying the relation between the parameters of the 
laser beams and the period and amplitude of the oscilla¬ 
tion. This is achieved by giving a set of formulae so that 
the relation can be understood analytically. The con¬ 


ditions that the frequency and/or the amplitude can be 
maximized or minimized have been predicted. To realize 
the mutual-transformation of the two components, not 
only the strength of the Raman coupling 7 is important, 
the energy difference between the two components is also 
important which can be tuned by 5 (this is similar to the 
case of a resonance). The frequency of oscillation will in¬ 
crease with 7 in general, while the amplitude depends on 
the interrelation of 7 and 5 . When (< 5 ,7) vary along a spe¬ 
cific straight line on the 6 — plane, the amplitude may 
increase, decrease, or even remain unchanged depending 
on the slope of the line. This is a common feature for 
the first case and for the second case with ^ = tt or 0. 
Whereas when 4 > tt and 0 , movable stripes will emerge 
on the ring. Experimental confirmation of the regularity 
unveiled in this paper and the predicted phenomena is 
expected. 
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Appendix 

When the initial state is ■i/'init = sm{(j)/2)(pq^ + 
cos((//2)(/3q4,, the associated up- and down-densities dur¬ 
ing the evolution are 


nt(6»,r) = ■;^{sin^((// 2 )[cos^(aq+/ 3 T)-h cos^( 2 p,+/ 3 )sin^(ag+/ 3 T)] 

ZTT 

-I- cos^{4>/2) sin^(2p,_^) sin^(a,_/ 3 T) 

— sin0sin(2pq_^) sin(ag_/3r)[cos(2/30 — 4,qj3T) cos{2pqj^p) sin(ag+/3T) 

-I- sin(2/30 - 49/3 t) cos(ag+/ 3 r)]}, (18) 

and 

n;(6»,r) = ^ {cos^ {(j) / 2)100^^ {aq-pr) + cos^ {2 pq_p)s\T?{aq-pT)] 

ZTT 

-h sin^((/)/2) s\T?{2pq+p) sin^(a,+/ 3 r) 

-I-sin0sin(2pq+^) sin(aq+/3r)[cos(2/30 — Agfir) cos{2pq-p) sm{aq-pT) 

-I-sin(2/30 — 4(7/3 t) cos(ag_/ 3 T)]}. (19) 
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